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The energetics of a Brownian motor and refrigerator driven by position dependent temperature, 
known as the Buttiker-Landauer motor and refrigerator, is investigated by extensive numerical 
simulations of the inertial Langevin equation. Our results are qualitatively different compared to 
previous works based on the overdamped Langevin equation and other phenomenological approaches. 
We find that the irreversible heat transfer via kinetic energy greatly reduces the efficiency /coefficient 
of performance of the motor/refrigerator. The motor and refrigerator work with maximum efficiency 
and coefficient of performance at optimal values of mass and friction coefficient but can never reach 
the Carnot efficiency /coefficient of performance. We also study the behavior of the motor in the 
linear response regime under condition of maximum power and find that efficiency at maximum 
power can never reach the efficiency of an endoreversible engine working at maximum power: the 
Curzon-Ahlborn efficiency. Under an analogous condition equivalent to that of maximum power, 
the coefficient of performance of the refrigerator cannot reach the coefficient of performance reached 
by a Carnot engine. Predictions of linear irreversible thermodynamics are in good agreement with 
numerical data. Finally, we investigate the role of different potential and temperature profiles to 
reduce the irreversible kinetic energy contribution, and increase the motor efficiency and refrigerator 
coefficient of performance. Our simulations show that optimizing the potential and temperature 
profile in order to reduce the irreversible heat transfer, diminishes the particle current as well 
leading only to a marginal enhancement of the system performance. 

PACS numbers: 05.10.Gg, 05.40.Jc, 05.70.Ln 



I. INTRODUCTION 

Technology has moved to an era where the creation 
and manipulation of devices at the nanoscale regime has 
become increasingly possible. At the nanoscale regime, 
the nanomachine is subject to thermal fluctuations which 
might seem to pose a hindrance for its functioning. How- 
ever, it is known that thermal noise can also play a ben- 
eficial role in the working of certain small scale devices. 
One example of such nanoscale devices is that of Brow- 
nian motors [l[ , which do useful work by rectifying ther- 
mal fluctuations. Since at equilibrium, the principle of 
detailed balance rules out any net flux, unidirectional 
motion can be achieved only when the system is driven 
away from equilibrium. 

The source of non-equilibrium could be an external sys- 
tem providing energy such as a time-dependent force or if 
the system is subject to a thermal gradient. Examples of 
the latter kind of Brownian motors are the autonomous 
thermal motors such as the Feynman-Smoluchowski (FS) 
motor [2||, where the system is simultaneously under the 
influence of two thermal reservoirs at different tempera- 
tures and the Buttiker-Landauer (BL) motor where 
the Brownian particle is alternately attached along the 
spatial coordinate, to thermal baths at different tem- 
peratures. Thermodynamic quantities such as the work 
output and efficiency, which characterize how effectively 
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such microscopic engines perform their tasks, have been 
a widely studied topic in recent years 

The energetics and transport properties of the BL mo- 
tor have been studied by many authors. Landauer 0], 
first proposed the idea of the BL motor when he showed 
the physical significance of nonuniform temperature in 
changing the relative stability of otherwise locally sta- 
ble states, a phenomenon he termed as the "blow-torch" 
effect. Such a position dependent temperature can also 
induce a net transport of Brownian particles in a peri- 
odic potential [1, [j, @| . Due of the asymmetry created 
by the nonuniform temperature, the Brownian particle 
can more easily climb the potential barrier when it is in 
a hot region than when it is in the cold region. As a 
result, in the steady state we find a net flux of Brownian 
particles alongwith a heat transfer in the same direction 
as the temperature gradient. In the presence of a small 
external load it can also do some work as a motor or heat 
engine. 

While a thermal gradient can create a net movement of 
Brownian particles, Onsager symmetry tells us that there 
will be a reciprocal process of the BL motor, namely the 
BL refrigerator. In presence of uniform temperature and 
an external force, there will be a heat flow from one ther- 
mal bath to the other, leading to cooling in one bath. In 
presence of a temperature difference, and an external load 
of magnitude greater than the stall load, the Brownian 
particles can transfer heat from the cold to the hot bath, 
like a heat pump @, H, [1[ . 

Stochastic energetics formulated by Sekimoto flOj ] . 
for mesoscale systems dominated by thermal fluctu- 
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ations, has been applied to investigate the efficiency 
of the BL motor/refrigerator. This led to a de- 
bate regarding the maximum attainable thermodynamic 
efficiency/coefficient of performance of the BL mo- 
tor/refrigerator. It was shown, using overdamped models 
that the BL motor can reach Carnot efficiency and the 
BL refrigerator can attain the corresponding Carnot co- 
efficient of performance 0, H, [Hj]. However, Derenyi- 
Astumian and Hondou-Sekimoto [IH argued that 
Carnot efficiency is unattainable due to the irreversible 
heat flow via kinetic energy. They made phenomeno- 
logical predictions regarding the failure of overdamped 
models in predicting this heat transfer. In a recent 
work [3], we confirmed Sekimoto's stochastic energetics 
for non-uniform temperature as well as the predictions 
of Derenyi-Astumian [l2j and Hondou-Sekimoto [llj], us- 
ing first principles molecular dynamics simulation and 
numerical solution of the Langevin equation, taking in- 
ertia into account. We also found that the irreversible 
heat transfer via the kinetic energy greatly reduces the 
efficiency of the motor. 

The kinetic energy contribution will also affect the re- 
frigerator performance in presence of non-uniform tem- 
perature as the irreversible heat transfer from the hot 
to the cold reservoir would reduce the amount of heat 
extracted from the cold bath thus diminishing the coef- 
ficient of performance. 

Various authors have studied the thermodynamics of 
the BL motor and refrigerator using overdamped mod- 
els @, [|| and other phenomenological approaches [§| . In 
the present work, we will discuss the energetics of the 
system by extensive numerical simulation of the inertial 
Langevin equation also compare the results in a few cases 
with first principles molecular dynamics simulation. We 
will also investigate the effect of different potential shapes 
and temperature profiles on the kinetic energy contribu- 
tion and whether the performance of the BL motor and 
refrigerator can be enhanced using this approach. 

Since a real motor must deliver a finite amount of work 
in a finite time, it is interesting from a practical viewpoint 
to also analyze its efficiency under the condition of max- 
imum power. In a recent work, Van den Broeck [TEj has 
shown that in the linear response regime, the efficiency 
at maximum power of a Brownian motor can attain the 
efficiency of an endoreversible engine when it operates at 
maximum power, namely the Curzon-Ahlborn efficiency 
if the condition of "tight coupling" is satisfied [16j]. The 
efficiency at maximum power of the BL motor has also 
been studied using overdamped models 0] and the au- 
thors showed that it is close to the Curzon-Ahlborn ef- 
ficiency. Corresponding to the efficiency at maximum 
power, some authors have studied the equivalent quan- 
tity for the refrigerator using linear irreversible thermo- 
dynamics [l7j . Since the overdamped model fails to pre- 
dict the correct heat transfer, we will examine the linear 
response behavior of the system by taking the inertia of 
the Brownian particle into account. 

The paper has been organized as follows: In the next 
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FIG. 1: (Color online) Two rectangular reservoirs filled with 
gas particles at temperatures Ti and T2 are alternately con- 
nected. Gas particles [red (dark gray) and blue (light gray) 
circles] are confined in the cell and only Brownian particles 
(large black circle) can move through the walls. Brown- 
ian particles are subjected to a piece-wise linear potential as 
shown in the bottom. 



section we introduce the basic model, as studied in our 
previous work [lij], and the methods used to investi- 
gate the energetics of the BL motor and refrigerator. In 
sections Hill and IIV| we discuss the efficiency and per- 
formance of the BL motor and refrigerator respectively 
based on the basic model. Section [V] discusses the ef- 
ficiency and coefficient of performance in the linear re- 
sponse regime. In section IVTI we study the effect of dif- 
ferent potential and temperature profiles on the motor 
efficiency and refrigerator coefficient of performance. Fi- 
nally, we end with a conclusion. 



II. BASIC MODEL AND METHODS 

We consider a chain of two-dimensional cells attached 
to each other along the x direction (see Fig. []}. Each 
cell is L/2 long and filled with N gas particles of mass 
to. The cells are thermally isolated from each other and 
at thermal equilibrium with a temperature independent 
from the neighboring cells. Brownian particles of mass M 
are placed in the cells and, unlike the gas particles, can 
move through the cell walls. They are also subjected to 
a periodic piece-wise linear potential U(x) with a period 
L in the x direction and constant in the y direction. For 
simplicity, we use a piecewise-linear potential: 



U(x) = 



'-^x for-f<a;<0, 
^■x for < x < f, 



(1) 



where Uo is the potential height. In addition to the peri- 
odic potential, a constant external force F is exerted on 
the Brownian particles. Instead of infinitely long chains 
we consider only two cells with a periodic boundary con- 
dition: cell 1 with temperature T\ for —L/2 < x < and 
cell 2 with T 2 for < x < L/2. 
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Apart from long-time fluid dynamical effects in two 
dimensions, the motion of the Brownian particle in the 
x direction can be investigated by the one-dimensional 
Langevin equation: 



(2) 



Mb = -~/(x)v - U'{x) +F+ y/2-y{x)T(x)£(t), 



where x and v are the position and velocity of the Brow- 
nian particle and £(t) is a standard Gaussian white noise: 

«(<)>= o, mm=s(t-s). (3) 

The piecewise constant temperature is given by: 

, , [t x for -k < x < 0, 
[T 2 for < x < y- 

Temperature is measured in energy unit (fee — !)• Here, 
and later on, an overdot refers to derivative taken with 
respect to time and a prime means derivative taken with 
respect to space. Assuming that Einstein's relation, 
namely, D{x) — T{x)/ r y{x) holds locally, the friction co- 
efficient also depends on the position in the same way as 
the temperature through the relation 



7(x) = pa B y/2'KT{x)m, 



(•5) 



where, p is the number density of gas particles in a cell. 
The above analytical expression for j(x) holds for an 
ideal gas 

When the relaxation time r = M/7 is much smaller 
than a mechanical time scale to — yj ML 2 /Uq, the iner- 
tia of the Brownian particle may be neglected, allowing 
the use of the overdamped Langevin equation. However, 
we confirmed in our previous work [I4j that the over- 
damped model fails in predicting the heat transfer when 
temperature is inhomogeneous. Therefore, we use the 
inertial Langevin equation @ and compare the results 
with molecular dynamics simulation. 

In our MD simulation, both gas and Brownian particles 
are hard disks of diameter a and <tb, respectively. A cell 
of size 250 x 400 (L — 500) contains 1000 gas particles so 
that p = 0.01, low enough to use the Langevin equation. 
For convenience, we assume m = 1, a = 1 and Uq = 1. 
See Ref. [l4[ for the detailed algorithm of the simulation. 

In order to investigate the thermodynamic properties, 
we use stochastic energetics formulated by Sekimoto [lfj, 
[T9L |2C| . The heat flux from the gas particles in the i-th 
cell to the Brownian particles is defined as, 
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where, (• • • ),; indicates ensemble average taken while the 
Brownian particles are located in the i-th cell. Using the 
Langevin equation ([2|) , we obtain the total heat flux as a 
sum of three terms: 



M d 
Ydt 

Qf E + Q 



(x 2 ). + {U'(x)x) i -F{x) i 

fi (7) 



where the first two terms on the rhs are the kinetic en- 
ergy and potential energy contribution to the heat flux, 
respectively, and the last term is the Joule heat. Work 
done on a Brownian particle by the external force F in 
each cell is given by, 



Wi = Fix) 



(8) 



and the total work done on the Brownian particle is given 
by, 



W = F{(x) 1 + (x) 2 ) = F(x) 



(9) 



where, = (1)2 — (%)/2, since width of both cells are 
the same. 

In the steady state, the net energy flux to the Brown- 
ian particles will vanish, and thus the energy gained by 
the Brownian particles in cell 1 will be cancelled by the 
energy loss in cell 2. Therefore, heat flux from cell 1 to 
cell 2 via the Brownian particles is defined as 



Qi^ 2 = Qi + W x = Ql 
= -Q 2 -W 2 = - 



KE 



■Qi 



PE 



(10) 



III. EFFICIENCY OF THE MOTOR 



When two cells have different temperatures (Ti > T 2 ), 
the Brownian particles in the higher temperature cell can 
reach the higher-potential energy region more often than 
those in the low-temperature cell. Hence, the Brownian 
particles tend to move from the hot to the cold cell over 
the potential barrier, resulting in a net current in the 
positive direction even in the absence of an external force 
(F = 0). In the presence of an external load (F < 0), the 
Brownian particles can do work against it behaving as 
a motor. In this section, we discuss the thermodynamic 
efficiency of the BL motor, defined as 



V 



-w 



Qi 



(11) 



In Fig. [51 we plot the motor velocity, rate of work out- 
put or the power —W, various components of the heat 
flow out of cell 1 at temperature Ti, and efficiency as a 
function of the external load. As the external load im- 
pedes the motion of the motor, the velocity decreases and 
vanishes at the stall load F sta a- When F = 0, the aver- 
age velocity of the Brownian particles is not zero, but the 
work is zero. At the stall load, since the motor velocity 
goes to zero, the work output again vanishes. At a certain 
value of F, between and F sta ii, the work output reaches 
a maximum i.e., the motor delivers maximum power. 

Now we come to the heat transfer. When the Brow- 
nian particle crosses one period to the right it extracts 
potential energy Uq from the higher-temperature cell at 
the rate (2Uq/ L){x)i and dissipates the same amount of 
potential energy to the cold cell. This heat flow via po- 
tential decreases with the load as the velocity decreases 
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FIG. 2: (a) Average velocity (v), (b) power -W, (c) Qf E 
and Qi , and (d) efficiency of the motor as a function of 
\F\L/2Uq. In (a), (b), (c) and, (d) the thin horizontal solid 
line represents y — 0. The parameter values are Tj — 0.7, 
T2 — 0.3, o\b = 5 and M/m — 5. The Carnot efficiency 
Vc = 0.57. 



and ultimately vanishes at F sta ii ■ At the quasistatic limit 
of zero motor velocity, this heat transfer via potential en- 
ergy is reversible. 

On the other hand, Fig. [2] shows that the heat flow via 
kinetic energy varies slowly with the load and is positive 
even in the quasistatic limit, showing that a net heat 
flows from the hot to the cold reservoir even when the 
motor does not deliver any work. As explained in our 
previous work [Ti| , when the Brownian particle crosses a 
temperature boundary from the hot to the cold side it dis- 
sipates heat due its kinetic energy, to the cold bath. This 



kinetic energy contribution will be proportional to the 
temperature difference between the hot and cold baths. 
Even when the velocity of the Brownian particles is zero, 
they are still crossing the temperature boundary back 
and forth. Therefore, this heat continues to flow at the 
stall load. Since the heat flows always from the hot cell 
to the cold cell, independent of the direction of the Brow- 
nian particle's movement, it is irreversible heat. 

The presence of this entropy producing irreversible 
heat carried by the kinetic energy of the Brownian parti- 
cle, makes it impossible for the BL motor to reach Carnot 
efficiency r\ c = 1 — T1/T2. In fact, Fig. [3] shows that the 
efficiency is a magnitude order lower than r\ c and dimin- 
ishes to zero at the quasistatic limit (at the stall load). 

Figure [3] shows the power, heat flow via kinetic en- 
ergy, and motor efficiency as a function of the tempera- 
ture difference between the two cells. Overdamped mod- 
els [7j , which ignore the kinetic energy contribution, pre- 
dict that the efficiency increases with AT and then sat- 
urates. However, when we take inertia into account its a 
different story. We find that at smaller AT the work out- 
put is low leading to a low efficiency while at large AT, 
the heat flow via kinetic energy is large, again reducing 
the efficiency. Figure |3J also shows that the kinetic energy 
is proportional to the temperature difference between the 
cells, since greater is the temperature difference between 
the cells, greater is the difference between the average ki- 
netic energy of the Brownian particle in the hot and cold 
cells. There is an optimum AT at which the efficiency is 
maximum but it is far lower than the Carnot efficiency. 

We had shown in our previous paper [lj| that Q KE 
diverges as M~ x / 2 in the overdamped limit (M — > 0). 
This can be understood as follows: the Brownian particle 
thermalizes over a length scale £ t h — VthT = vTM/f, 
where v t h — ^JT/M is the thermal velocity. Within this 
thermal length scale, the average kinetic energy of the 
Brownian particle deviates from the bath temperature. 
If, we are in the overdamped regime and the width of 
the cells is much larger than the thermal length scales 
then, away from the border region betwen the hot and 
cold baths the kinetic energy of the Brownian particle 
is adapted to the local bath temperature [l3j]. However, 
in the narrow transition region between the two baths 
the average kinetic energy of the particle will gradually 
decrease from Ti/2 inside the hot bath to T2/2 inside 
the cold bath (see the kinetic energy profile in Fig. 3 of 
Ref. 0). 

The kinetic energy contribution, as Fig.[3]shows, is pro- 
portional to the temperature difference between the cells. 
As the temperature is non-uniform, the kinetic energy 
contribution will consequently be proportional to the ef- 
fective temperature gradient: | T\ — T2 | jlth oc M ~ x / 2 . 
As we approach the overdamped limit, the length of the 
transition region decreases and the average kinetic energy 
profile of the Brownian particle in this transition region 
becomes more and more steep. In the overdamped limit 
[M — > 0), the change in the average kinetic energy will 
be sudden as the transition region vanishes and hence the 
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FIG. 3: (a) Power, (b) heat flux due to kinetic energy, and (c) 
efficiency as a function of AT — Ti — T%, where T2 is kept at 
0.1. The parameter values are M/m = 5.0, \F\L/2U = 0.5 
and <jb = 5.0. In (a), (b) and, (c) the thin horizontal dashed 
line represents y = 0. The Carnot efficiency varies from 0.67 
at AT = 0.2 to 0.96 at AT = 2.8. 



effective thermal gradient tends to infinity and the irre- 
versible heat transfer diverges as M~ x l 2 . Since this irre- 
versible heat dominates the heat flow via potential, the 
efficiency diminishes to zero as VaJ, in the overdamped 
limit. This was also predicted by Derenvi-Astumian[12j| . 

On the other hand, in the underdamped limit (M — ► 
00), the Brownian particle doesn't thermalize to the local 
environment before entering the next heat bath [2lT | and 
the motor fails, diminishing the power and efficiency to 
zero. The efficiency will be maximum at an optimal value 
of M. Figure [U clearly confirms these predictions. We 
also find good agreement between the inertial Langevin 
equation and molecular dynamics simulation, suggesting 
the validity of Sekimoto's stochastic energetics for non- 
uniform temperature. 

We now study the effect of friction coefficient on the 
energetics. In Fig. [5] we plot the power, heat flow via 
kinetic energy and efficiency as a function of the friction 
coefficient. In the underdamped regime (7(1) <C 1), the 
kinetic energy of the Brownian particle does not equili- 
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FIG. 4: (a) Power, (b) heat flux due to kinetic energy, and (c) 
efficiency as a function of the mass ratio M/m. Symbols and 
lines correspond to MD simulation and the inertial Langevin 
equation respectively. The dotted line in (b) corresponds to 
a phenomenological fit (oc M _1//2 ) to the kinetic energy con- 
tribution. The dashed line in (c) represents a fit (oc i/M) to 
the efficiency in the overdamped limit [M -t 0), In (a), (b) 
and, (c) the thin horizontal solid line represents y = 0. The 
parameter values are Ti = 0.7, T 2 = 0.3, \F\L/2U = 0.3 and 
(jb = 5.0. The Carnot efficiency r\ c — 0.57. 



brate to the bath temperature and hence the asymmetry 
due to the inhomogeneous temperature becomes negligi- 
ble. As a result, a positive current cannot be sustained 
against the load, and the motor fails to deliver work. 
Similarly, in the overdamped limit (7(1) — > 00) the large 
friction coefficient makes the Brownian particle immobile 
and the work output again goes to zero. 

However, the irreversible heat flow via kinetic energy 
shows a different behavior as the friction coefficient in- 
creases. As the medium becomes more damped, the 
Brownian particle tends to adapt to the local bath tem- 
perature so that its average kinetic energy in each cell 
approaches Tj/2 (i = 1,2). Since, the kinetic energy con- 
tribution Q™ 2 > i s proportional to the difference between 
the average kinetic energy of the Brownian particle in the 
two cells, it increases with the damping and saturates 
in the overdamped limit, confirming Hondou and Seki- 
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FIG. 5: (a) Power, (b) heat flux due to kinetic energy and, 
(c) efficiency as a function of 71. The parameter values are 
Ti = 0.7, T 2 = 0.3, \F\L/2U = 0.2 and M/m = 5.0. In (a), 
(b) and, (c) the thin horizontal dashed line represents y — 0. 
The friction coefficient is varied by increasing the diameter of 
the Brownian particle. The Carnot efficiency 7/ c = 0.57. The 
top x — axis in (a) shows the dimensionless friction coefficient 
'yiL/y/MUo, as an indicator of the friction regime. Typically, 
jiL/y/MUo > 10 corresponds to the overdamped regime and 
"fiL/y/MUa < 1 corresponds to the underdamped regime. 



moto's conjecture that the kinetic energy contribution is 
independent of the friction coefficient in the overdamped 
regime. Finally, Figure [5] shows that there is an optimum 
friction coefficient at which the efficiency is maximum. 



IV. PERFORMANCE OF REFRIGERATOR 

We now discuss the the energetics of the BL refriger- 
ator. When the external load exceeds the stall load, the 
motor will move backward in the same direction as F. In 
this situation, some authors 0, @, 0] have claimed that 
the Brownian particle will transfer heat against the tem- 
perature gradient, from the cold to the hot bath, func- 
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<3™)i an d (c) the coefficient of performance e as a function 
of F/Fstall, F s taii being the stall load. In (a), (b) and, (c) the 
thin horizontal dashed line represents y — 0. The parameter 
values are AT = Ti - T 2 = 0.01, T 2 = 0.5, M/m = 5.0 and 
ob = 5.0. The Carnot coefficient of performance e c = 50. 



tioning as a refrigerator, and at the stalled state the co- 
efficient of performance e, defined as, 



Q2 



(12) 



will attain the Carnot coefficient of performance e c = 
T 2 /(T 1 -T 2 ). 

In Fig. O we plot the heat flow out of the hot and cold 
cells as a function of the external load. In contradiction 
to the claim made in Refs. @, B @, we find that the 
stalled state, at which the current changes its direction 
is not the transition point from the motor to the refrig- 
erator. At the stalled state, the heat flow via potential, 
which follows the direction of the particle flux (v), goes 
to zero. However, there is still a net heat transfer from 
the hot to the cold cell due to the kinetic energy contri- 
bution. In fact, the heat flow via kinetic energy never 
vanishes and always flows from the hot to the cold reser- 
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bottom the heat flows are 

QT, Q2, Qf E and Qi) and, (c) the 
coefficient of performance e as a function of AT for M/m = 
5.0, ob = 5.0, \F\L/2U = 0.5. In (a), (b) and, (c) the thin 
horizontal dashed line represents y — 0. Temperature of the 
cold bath is kept constant at T2 = 0.5 and Ti is varied. e c 
varies from 00 at AT = to 20 at AT = 0.025. 



FIG. 8: (a) W, (b) various components of Q2 (from top to 
bottom the heat flows are 

Q™ Q2, QP and Qi) and, (c) 
the coefficient of performance e as a function of M/m for 
AT = 0.01, a B = 5.0, \F\L/2U = 0.5, Ti = 0.51, T 2 = 0.5. 
In (a), (b) and, (c) the thin horizontal dashed line represents 
y — 0. The Carnot coefficient of performance is e c = 50. 



voir regardless of the direction of particle current and 
magnitude of the external load. 

Figure [5] shows that starting from the stall load there 
is a region where the model works neither as a motor nor 
as a refrigerator i.e., heat flows from the hot to the cold 
cell and work is done on the Brownian particle. Refriger- 
ation starts only when the heat flow via potential exceeds 
the sum of joule heat dissipated to the cold bath and the 
heat transfer via kinetic energy. As the load is further en- 
hanced, the joule heat (varying as |T| 2 ), which originates 
from the power supply W, eventually dominates the heat 
flow via potential (varying as \F\) and the cooling stops. 
There is an optimal load at which the heat flow out of 
the cold cell and the coefficient of performance attain a 
maximum but e is far below e c , due to the heat flow via 
kinetic energy. 

In Fig. we show the power supply, various compo- 
nents of the heat transfer from the cold bath at temper- 



ature T2 and the coefficient of performance as a function 
of the temperature difference between the cells. We keep 
T2 constant and vary T±. As the temperature difference 
between the baths increases, the net negative current in 
the direction of F will decrease as the increasing temper- 
ature of the hot bath will try to drive the Brownian parti- 
cle in the forward direction. This decreases the heat flow 
via potential, which is proportional to the current, and 
hence the refrigerator's ability to transfer heat against 
the thermal gradient. The power supply and joule heat 
also decrease with AT. However, as AT <C 1, this de- 
crease is marginal. On the other hand, the irreversible 
heat dumped to the cold cell via kinetic energy, increases 
significantly with AT and greatly reduces the cooling. 
In fact, beyond a small temperature difference AT <C 1, 
refrigeration fails, as Fig.[7]shows. The coefficient of per- 
formance is monotonically reduced with increasing tem- 
perature difference between the cells. 

Figure [5] shows the power input, various components of 
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FIG. 9: (a) W, (b) various components of Q2 (from top to 
bottom the heat flows are , Q2, Q2 an d Q2) an d, (c) the 
coefficient of performance e as a function of 71 for AT = 0.01, 
M/m = 5.0, \F\L/2U = 0.5, 7\ = 0.51, T 2 = 0.5. Carnot co- 
efficient of performance e c = 50. In (a), (b) and, (c) the thin 
horizontal dashed line represents y — 0. The friction coeffi- 
cient 71 is varied by increasing the diameter of the Brownian 
particle. The top a;— axis in (a) shows the dimensionless fric- 
tion coefficient jiL/y/MUo, as an indicator of overdamping. 



the heat flow and the coefficient of performance as a func- 
tion of the mass of the Brownian particle. The larger the 
mass of the Brownian particle the larger is the intrawell 
relaxation time thus leading to a low particle current in 
the negative direction. Hence the heat flow via poten- 
tial as well as the power input and joule heat decrease 
with the inertia of the Brownian particle. As expected, 
Fig. [8] shows that the irreversible heat via kinetic energy 
diverges in the overdamped limit. There is another com- 
ponent of the kinetic energy contribution, which goes to 
zero in the overdamped limit but becomes significant for 
large inertia. In the overdamped regime, the change in 
potential energy is immediately dissipated to the bath as 
the Brownian particle goes down the slope. However, at 
finite mass, the potential energy change is transported to 



the next cell as an "unthermalized" kinetic energy con- 
tribution, where it is dissipated as heat. This "unther- 
malized" kinetic energy contribution increases with the 
mass of the Brownian particle Hence, the total heat 
flow and the coefficient of performance go to zero in both 
the underdamped and overdamped limits. There is an 
optimal mass of the Brownian particle at which Q2 and 
e are maximized. But, e is far below the Carnot limit. 

In Fig [5] we plot W, various components of Q2 and 
e as a function of the friction coefficient. Since our nu- 
merical algorithm for solving the Langevin equation fails 
when 7(x)L/VA/TJo -> 0, we could not reach the ex- 
treme underdamped limit j(x) — > 0. In the underdamped 
regime, the "unthermalized" component of the kinetic en- 
ergy is large, decreasing the total heat flow Q2- In the 
overdamped limit, the kinetic energy contribution due to 
the temperature difference between the cells, saturates 
at a certain value. However, the large friction coefficient 
makes the Brownian particle almost immobile thus reduc- 
ing the particle current and consequently the heat flow 
via potential. The work input and joule heat decrease as 
well. Hence, the coefficient of performance again dimin- 
ishes to zero in both the overdamped and underdamped 
limits. There is an optimum friction coefficient at which 
Q2 and e are maximum. However, e is far below e c . 

In our previous work flij , we had shown that when 
two cells have the same temperature and the Brownian 
particles are driven by an external force F > 0, we ob- 
serve a net cooling in cell 1 at the expense of heating of 
cell 2. In Fig. [TU1 we show components of the heat flow 
from cell 1 to the Brownian particle and the coefficient of 
performance e = Q\/W as a function of F when T\ = T?. 
We find that there is an optimal F at which Qi is max- 
imized. When temperature is uniform we do not expect 
any kinetic energy contribution. Hence it follows that e 
will reach e c , which goes to 00 when T\ =T%. Figure [TUb 
shows that e tends to infinity as F — > which means that 
when temperature is uniform, a vanishingly small exter- 
nal force (or a power supply W which goes to zero) is 
sufficient to generate a heat flow from cell 1 to cell 2. We 
also find good agreement between the inertial Langevin 
equation and the molecular dynamics simulation. 



V. LINEAR IRREVERSIBLE 
THERMODYNAMICS 

In this section we discuss the efficiency and coefficient 
of performance of the BL system in the linear response 
regime. The BL motor and refrigerator are a result of 
cross effect between the external force F and the tem- 
perature difference AT = T\ — T%. If they are small, 
there exists a linear relationship between the thermody- 
namic fluxes ((v) and Qi_>2) and thermodynamic forces 
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FIG. 10: (a) Work input, (b) various components of Qi of the 
refrigerator as a function of FL/2Uo, (from top to bottom, the 
heat flows are Qi E , Qi, Q™ and Q\), and (c) e as a func- 
tion of FL/2Uo, obtained from the inertial Langevin equa- 
tion (lines) and molecular dynamics simulation (symbols). In 
(a), (b) and, (c) the thin horizontal dashed line represents 
y = 0. The parameter values are T = 0.5, M/m = 20.0 and 
a B = 8.0. 



(F/T and AT/T 2 ): 



Qi 



F 
T 



= L 



21 



T 



AT 
'■ J^2 ' 

AT 

' 



(13) 



^22' 



where Ly's — 1,2) are the Onsager transport co- 
efficients, Ti = T + AT/2 and T 2 = T - AT/2. While 
investigating the Biittiker-Landauer motor and refrigera- 
tor, we had shown in our previous work [l4l |. the validity 
of linear irreversible theory and numerically confirmed 
Onsager symmetry relationship, namely L 12 = T 2 i, for 
spatially inhomogeneous temperature. 
The efficiency of the motor is given by, 



-W 



—F ( Lu — 



J\2- 



AT 



yi^2 ^2lTp + L22^j?r 

Recently, Van den Broeck [z| investigated the efficiency 
of Brownian motors in the linear response regime and 



concluded that, in principle, Carnot efficiency can be at- 
tained. He argued that for small F, the temperature dif- 
ference AT drives the motor in the forward direction and 
simultaneously there is a heat transfer from the hot to 
the cold reservoir. As F increases, we ultimately reach 
the quasistatic condition of zero motor velocity at the 
stall load F sta u- From Eq. (fTB"]) . we find that the stall 
load is given by, 



F. 



stall 



_Lu / AT 
Ln { T 



(15) 



Beyond the stall load, the motor reverses its direction, 
transferring heat against the temperature gradient as a 
heat pump. In an ideal system, the velocity and heat 
reverse directions at the same value of F sta ii so that (v) 
and Qi_*2 vanish for nonzero F and AT. Substituting 
the value of F sta ii for F in Eq. (fT4)) . we find that this 
happens only when the condition of zero entropy pro- 
duction: TnT 2 2 — L12L21 holds, and Onsager symme- 
try relationship is valid. This condition is referred to as 
"tight coupling" condition and when it is satisfied, we 
find that efficiency at the quasistatic limit is equal to the 
Carnot efficiency: AT/T. 

Carnot efficiency, however, can only be reached in the 
quasistatic condition where the motor velocity is zero, 
and hence it does not deliver any work. In practice 
a motor, must perform finite work in a finite time. A 
more practical measure of the motor efficiency is the ef- 
ficiency when it delivers maximum work. Curzon and 
Ahlborn [Til ] claimed that the efficiency at maximum 
power 77* is bounded by the Curzon- Ahlborn efficiency, 
Vmax = 1 ~ ^/T 2 /Ti, where T\ > T 2 . The claim 
was found to be valid for thermal motors Asfaw 
and Bekele Q studied the efficiency at maximum power 
(77*) of the BL motor and found that i]* approaches 
Vmax as the temperature difference between the reser- 
voirs T\ — T 2 = AT — > 0. However, their conclusions 
are not reliable since they used the overdamped model 
and also a different definition of work from the one used 
in the Carnot efficiency. In order to evaluate efficiency 
at maximum power, we must use the proper thermody- 
namic definition of work (Eq. @) and take into account 
the inertia of the Brownian particle into account. 

Maximum power will be attained at a value of the load 
F* given by, 



= _Lu / AT 
Tn V 2T 



(16) 



Substituting the value of F* in Eq. (|T4|) , we find that the 
efficiency at maximum power is given by, 



where, 



AT q 2 
~2T2-q 2 



L\2 

V TnT 22 



(17) 



(18) 
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If the necessary condition of zero entropy production is 
satisfied, then q = 1 and to the lowest order of AT, the 
efficiency at maximum power approaches the efficiency 
of an endoreversible heat engine at maximum power viz., 
the Curzon-Ahlborn efficiency: rj* nax = 1 — \/To/T-\ ~ 
AT/2T [l||. However, in our previous work [14| we 
had shown that the product of the diagonal coefficients 
T11T22 diverges as (M/m)~ 1//2 , reflecting the divergence 
of the kinetic energy contribution. The off-diagonal co- 
efficients however, did not show this singular behavior. 
Consequently, we find L11L22 £12-^21 for all values 
of M/m indicating large entropy production and very 
low efficiency. Hence the BL motor can never attain the 
Carnot efficiency or the Curzon-Ahlborn efficiency. 

The coefficient of performance of the refrigerator is 
given by, 



e = 



-Q 



l->2 



^21" 



t AT 
-^22 j^ - 



W 



F in 



t AT 



(19) 



By similar arguments as in the case of the motor, it is 
easy to show that the coefficient of performance reaches 
the Carnot coefficient of performance TJ AT in the qua- 
sistatic limit, only if the necessary condition -L11L22 = 
Z/12-L21 holds. 

In the quasistatic limit, the refrigerator does not trans- 
fer any heat and hence is practically useless. Since 
Carnot coefficient of performance is only possible in the 
quasistatic limit, it is not an effective measure of the re- 
frigerator performance from a practical point of view. A 
refrigerator must transfer a finite amount of heat in a 
finite time. Similar to the efficiency at maximum power 
for the motor we must find the coefficient of performance 
of the refrigerator, under conditions of finite heat trans- 
fer from the cold to the hot bath. Following Van den 
Broeck's work fl5| , regarding the efficiency at maximum 
power for the motor, some authors have investigated the 
coefficient of performance of the refrigerator at maximum 
Qi^ 2 (AT/T) [13. This quantity is maximum at a tem- 
perature difference AT* , given by 



AT* 



L21TT 
£22 2 



(20) 



Substituting the value of AT* in place of AT in Eq.(fl9]). 
we find that up to the lowest order of AT, the coefficient 
of performance at maximum Qi_>2(AT/T), is given by 



T 



q- 



2AT 



(21) 



Analogous to the efficiency at maximum power, e* is 
equal to a factor that depends on q only, times half the 
Carnot coefficient of performance. If the necessary condi- 
tion for attaining the Carnot coefficient of performance, 
namely L11T22 = £12^21 holds, then q = 1 and e* attains 
a maximum e* max given by, 



T 
2AT 



(22) 




FIG. 11: rf (symbols), rjg (dashed line) and the Curzon- 
Ahlborn efficiency i;^, (solid line) as a function of AT = 
Ti — T2, where the temperature of the cold cell is kept at 
T2 = 0.5. ?7*naa; has been scaled (i)^ OI x 0.075) for easy 
visualization of all quantities on the same plot. The inset 
shows the detail of the efficiency at maximum power in the 
linear response regime. Parameter values are M/m = 5.0 and 
ob = 0.5. 



In terms of the load, we can write this expression as, 

£21 1 



Lii2LF| 



(23) 



However, for the BL refrigerator, q = 1 will never be 
attainable due to the kinetic energy contribution [lij and 
consequently e* will be less than ££, os - 

These predictions of linear irreversible thermodynam- 
ics concerning rf and e* have not been verified when 
temperature is spatially inhomogeneous. One of our ob- 
jectives is to investigate the validity of linear irreversible 
theory via numerical simulation of the inertial Langevin 
equation. 

From Fig. [2] we find that the power —W reaches a 
maximum at an optimum value of the load F*. We then 
evaluate the efficiency at F* , which is the efficiency at 
maximum power, as a function of AT. 

In Fig. [TTJ we plot as a function of AT, the Curzon- 
Ahlborn efficiency rf obtained from numerical so- 
lution of the inertial Langevin equation and the efficiency 
at maximum power predicted by linear response theory, 
namely rf q . In order to determine rf q we have to find q, 
which is evaluated from the values of the Onsager coeffi- 
cients, numerically. As the kinetic energy contribution is 
proportional to the temperature difference, it diminishes 
rf to zero as AT — > 00, thus contradicting the results of 
Ref. 0, obtained using overdamped models. There is an 
optimum AT at which rf reaches a maximum but it is 
far below f^ax- I n ^he linear response regime, there is 
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FIG. 12: e* obtained from numerical simulation (symbols), 
e* q (dashed line) and e* ma:c (solid line), as a function of the 
external force. Parameter values are M/m = 5.0, as = 5.0 
and 7i = T 2 = 0.5. 
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FIG. 13: Heat flow via kinetic energy at the potential minima 
(dashed line) and maxima (solid line) as a function of the 
external load. Parameter values are Ti = 0.7, T2 = 0.3, 
M/m — 5 and as = 5.0. 



good agreement between rj* and rj* , as shown in the inset 
of Fig. [dp. 

We now investigate the performance of the refrigera- 
tor at maximum Q2-+1AT/T , numerically. The quantity 
Q2-+1 AT/T reaches a maximum at an optimum value 
of the temperature difference, namely AT*. Then, we 
evaluate the coefficient of performance e* at AT* as a 
function of |F|. Figure [T2l shows that e* is always lower 
than £* max due to the kinetic energy contribution. There 
is also good agreement between e* and the coefficient of 
performance predicted by linear response theory e*. 

VI. EFFICIENCY OPTIMIZATION 

Figure 2] ([5]) showed that 77 (e) reaches a maximum at 
optimum values of the mass and damping between the 
overdamped and underdamped limits. In this regime, 
the efficiency can be further enhanced by reducing the 
kinetic energy contribution. 

The major cause of the low efficiency and coefficient 
of performance in the BL motor and refrigerator respec- 
tively is the heat transfer due to kinetic energy across 
the temperature boundary. When a temperature bound- 
ary coincides with the potential minima, the situation is 
the worst since the Brownian particle crosses the tem- 
perature boundary many times, resulting in a large heat 
transfer via kinetic energy, as illustrated in Fig. 1131 In 
order to reduce the number of crossings, we must move 
the temperature boundary away from the potential min- 
ima. On the other hand, the other temperature boundary 
should remain close to the potential maxima to maintain 
the current of Brownian particles. 
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FIG. 14: (Color online) The temperature and potential pro- 
files corresponding to case I of motor and refrigerator. The 
dark circle represents the Brownian particle. 



Another way to minimize the number of crossings 
would be to keep the same temperature profile but mod- 
ify the potential shape so that there is a small potential 
barrier of height E < Uq at the location of the poten- 
tial minima, which in our original model coincides with 
the temperature boundary. Since the Brownian parti- 
cle is less likely to be found at a potential maxima, the 
number of crossings will decrease and consequently the 
kinetic energy contribution will become diminish. 

Based on the above ideas, we will study the optimiza- 
tion of the BL motor and refrigerator in two different 
cases, where we vary the potential (case /) and tempera- 
ture {case IT) profiles as compared to our original model. 



A. Motor 

In case I, we shift the location of the temperature 
boundary away from the potential minima. The poten- 
tial U{x) is the same as in our original model but the 
temperature profile T(x) is modified and is now given 
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FIG. 15: (a) Power (— W), (b) heat flux due to kinetic energy, 
and (c) efficiency of the motor as a function of 5/L for case I 
of motor obtained from the inertial Langevin equation. In (a) , 
(b) and, (c) the thin horizontal dashed line represents y = 0. 
Parameter values are Ti = 0.7, T 2 = 0.3, \F\L/2U = 0.1, 
M/m = 5 and ob = 5.0. The Carnot efficiency rj c = 0.57. 



by: 



T(x) 



To for 



L 



< x < 5 



12 iU1 ~~2 

Ti for 5 < x L 



(24) 



where, T\ > T2 and the parameter <5 specifies the location 
of the temperature boundary. 

At 5/L — 0, the cell boundary coincides with the po- 
tential maxima/minima and we go back to our previous 
model studied in section [TTJ At 5/L = —1/2 and 1/2, the 
temperature becomes uniform since the width of the cold 
(5/L = -1/2) /hot (5/L = 1/2) cell goes to zero and the 
system ceases to work as a motor. The motor can deliver 
work for 5/L in the range -1/2 < 5/L < 1/2. In Fig.[T5l 
we plot the work output, heat flow and efficiency as a 
function of the parameter 5/L. Figure [T5l clearly shows 
that as the cell boundary moves away from the potential 
minima, in either direction, the irreversible heat trans- 
fer via kinetic energy is reduced. When 5/L > 0, the 
presence of the cold region reduces the likelihood of the 
Brownian particle jumping over the barrier in the positive 




FIG. 16: (Color online) The temperature and potential pro- 
files corresponding to case II of motor and refrigerator. The 
dark circle represents the Brownian particle. 



direction, especially if (2Uo/L)5 3> T2, thus diminishing 
the current and the power —W, and consequently the ef- 
ficiency. On the other hand if 5 /L is close to the potential 
minima, the kinetic energy contribution becomes large, 
again reducing the efficiency. At an optimum 5/L > 0, 
the efficiency is maximized as shown in Fig. 1151 

When 5/L < 0, the hot region extends to the left 
of the potential minima at x = 0. This enhances the 
probability of the Brownian particle reaching the higher- 
potential-energy region and crossing the barrier in the 
negative direction. As a result, the net particle current 
in the positive direction and the power delivered by the 
Brownian particle diminish, thus reducing the efficiency. 
On the other hand, if 5/L is close to x = 0, the large 
heat transfer via kinetic energy reduces the efficiency. 
Figure [TBI shows that there is an optimum 5/L < 0, at 
which the efficiency is maximized. However, the maxi- 
mum efficiency is far below the Carnot limit. 

In case II, we keep the same temperature profile as 
in our original model [Eq. ([4])] but change the potential 
profile U(x), which is now given by, 



2Ug 
L-2A 



(x + A) for-^<x^-A 



U(x) 



f(* + A) 

-f(z-A) 

If^(z-A) 



for -A < x ^ 
for < x ^ A 
for A < x ^ y 



(25) 



Equation( [25]) shows that the potential has a barrier of 
height E < Uq at x = 0, coinciding with the position of 
the potential minima in our original model. 

In Fig.[T71 we plot —W, various components of Qi and 
X] as a function of the barrier height E at x = normal- 
ized by Uq. At E/Uo — 0, the temperature boundary 
coincides with a potential minima and the kinetic energy 
contribution is large thus leading to a low efficiency. As 
we increase E/Uq, the irreversible heat transfer via ki- 
netic energy decreases as the recrossing of the Brownian 
particle near the temperature boundary at x = is re- 
duced on account of a potential maxima there. However, 
the motor velocity and consequently the power is dimin- 
ished as well because the asymmetry around a poten- 
tial minima due to the inhomogeneous temperature de- 
creases. The efficiency is maximized for < E/Uq < 1, 
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FIG. 17: (a) Power, (b) heat flux due to kinetic energy, and 
(c) efficiency of the motor as a function of E/Uo for case 
II of motor. Solid lines correspond to data obtained from 
the inertial Langevin equation. Circles correspond to data 
obtained from MD simulation. In (a), (b) and, (c) the thin 
horizontal dashed line represents y — 0. Parameter values are 
Tj = 0.4, T 2 = 0.2, A/L = 0.1, \F\L/2U = 0.025, M/m = 5 
and cjb = 5.0. The Carnot efficiency rj c = 0.5. 



FIG. 18: (a) Work input W, (b) various components of Q2 
(from top to bottom, the heat flows are Q™! Q2, Q\ and 
Q™), an d (c) the coefficient of performance e of the refriger- 
ator as a function of S/L, corresponding to case I of refrig- 
erator obtained from the inertial Langevin equation. In (a), 
(b) and, (c) the thin horizontal dashed line represents y = 0. 
The parameter values are \F\L/2Uq = 0.2, M/m = 5.0 and 
ob = 5.0, Ti = 0.51 and T 2 = 0.5. The Carnot coefficient of 
performance e c = 50. 



as shown in Fig. 1171 though it remains far below the 
Carnot limit (rjc = 0.5). Figure [171 also shows qualita- 
tive agreement between the inertial Langevin equation 
and MD simulation. 



B. Refrigerator 

We now consider the performance of the BL refriger- 
ator corresponding to the same two cases, which were 
introduced in the previous subsection for the BL motor. 

In Fig. [TH1 we plot the work input, heat flow out of cell 
2 and the coefficient of performance as a function of the 
location of the temperature boundary 5/L correspond- 
ing to case I. Since AT <C 1, the thermal gradient has 
a negligible effect on the particle current. Instead, the 
magnitude of the negative current is determined mostly 
by the potential shape and the external force F. How- 



ever, the potential profile stays the same regardless of the 
position of the temperature boundary. Hence, the power 
input varies slowly with 5/L, as shown in Fig. 1181 

As expected, Fig. [15] shows that the kinetic energy con- 
tribution becomes smaller as the temperature boundary 
is shifted away from the potential minima. The joule heat 
dissipated to the cold cell vanishes at S/L = —1/2 since 
the width of the cold region goes to zero at that point. As 
5/L shifts to the right of x = —1/2, the joule heat starts 
to increase as the width of the cold cell becomes larger 
and reaches a maximum at 5/L = 1/2, when the cold cell 
occupies the entire period L of the potential. As the tem- 
perature boundary is shifted to the left of 5/L = 0, the 
heat flow via potential Q™, is reduced since the width 
of the cold cell becomes smaller and less potential en- 
ergy is extracted. When 5/L > 0, Qf E again diminishes 
since the cold cell now includes a portion of the negative 
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slope as well. Figure [19] shows that there is an optimum 
6/L < at which the total heat extracted from cell 2 and 
the coefficient of performance reach a maximum, showing 
that this model can enhance the cooling and refrigerator 
performance by appropriately changing the location of 
the temperature boundary. However, e is still far below 
the Carnot coefficient of performance. 

In case II, the small barrier at x = reduces the heat 
transfer via kinetic energy. This is clearly seen in Fig. [TUl 
which shows the work input, heat flow out of cell 2 and 
e as a function of the barrier height E at x = 0, normal- 
ized by Uq. Due to the small temperature difference, the 
current and consequently the work input and the joule 
heat vary slowly with E/Uq. On the other hand, the 
barrier of height E at x = diminishes the heat flow via 
potential because due to the negative slope, —E/A, an 
amount of heat E is dissipated each time the Brownian 
particle crosses the cold cell, thus reducing the magni- 
tude of potential energy extracted to Uq — E. Since, heat 
flow via potential decreases with E/Uq more rapidly as 
compared to the kinetic energy contribution, there is no 
optimum value of E/Uq and net heat extracted from cell 
2 decreases monotonically from a maximum at E/Uq = 0. 
The coefficient of performance also decreases with E/Uo 
and is far below e c . Hence, this case does not improve 
the performance of the refrigerator. 



VII. CONCLUSION 
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We studied the efficiency and coefficient of perfor- 
mance of the Buttiker-Landauer motor and refrigerator 
by numerical simulation of the inertial Langevin equa- 
tion and MD simulation. Our results show qualitatively 
different behavior of the heat transfer, efficiency and co- 
efficient of performance as a function of various parame- 
ters as compared to previous works based on overdamped 
models 0, H[ or other phenomenological approaches @ . 
From extensive numerical simulations, we find that the 
irreversible heat flow via kinetic energy greatly reduces 
the motor efficiency and refrigerator coefficient of perfor- 
mance. As long as there is a temperature difference be- 
tween the heat baths, this heat flow cannot be neglected. 
As a consequence, the transition point from motor to re- 
frigerator is not at the stalled state, but when the heat 
transfer via potential is greater than the sum of the joule 
heat and the heat flow via kinetic energy. This occurs 
only when the load exceeds the stall load. 

There is an optimal value of the mass and friction coef- 
ficient at which the efficiency /coefficient of performance 
reach a maximum. However, they are far below their 
respective Carnot limits. We also found that the effi- 
ciency at maximum power can never attain the Curzon- 
Ahlborn efficiency due to the irreversible heat flow via 
kinetic energy. Due to this irreversible heat, the coef- 
ficient of performance of the refrigerator is below the 
theoretical maximum predicted by linear response the- 
ory, under conditions which could be considered as the 



FIG. 19: (a) Work input W , (b) various components of Q2 
(from top to bottom, the heat flows are Q™! Q2, Qi and 
<32 E )i an d (c) the coefficient of performance e of the refriger- 
ator as a function of E/Uo, corresponding to case II of refrig- 
erator obtained from the inertial Langevin equation. In (a), 
(b) and, (c) the thin horizontal dashed line represents y = 0. 
The parameter values are \F\L/2Uo ~ 0.2, M/m = 5.0 and 
a B = 5.0, Ti = 0.51 and T 2 = 0.5. The Carnot coefficient of 
performance e c = 50. 



equivalent of maximum power. The predictions of lin- 
ear irreversible thermodynamics are in good agreement 
with numerical data. Changing potential shape or the 
temperature profile in order to reduce the heat flow via 
kinetic energy reduces the particle current as well, lead- 
ing only to a small enhancement of the motor efficiency 
and refrigerator coefficient of performance. 

Since reduction of heat transfer via kinetic energy re- 
sults in diminishing the particle current thus making the 
motor less efficient, the next question to address is how 
one can reduce the irreversible heat without decreasing 
the particle current in order to enhance the motor effi- 
ciency and refrigerator coefficient of performance. At a 
more fundamental level, it would be interesting to obtain 
theoretically the maximum efficiency at maximum power 
and the equivalent maximum for the coefficient of per- 
formance. Answering such questions would also enhance 
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our understanding of non-equilibrium thermodynamics 
at the small scale. 
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